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SPACE STATION/ SHUTTLE ORB ITER 
DYNAMICS DURING DOCKING 

N. G. Fitz-Coy* and J. E. Cochran, Jr. + 

Auburn University, AL 36849-3501 

Mathematical models of a reference Space Station 
configuration ("Power Tower") and a Space Shuttle 
Orbiter are developed and used to study the dynamic 
behavior of the Space Station/Orbiter system just 
prior to and subsequent to an impulsive docking of 
the two spacecraft. The physical model of the Space 
Station is a collection of rigid and flexible bodies. 
The orbiter is modeled as a rigid body. An algorhthm 
developed for use in digitally simulating the dynamics 
of the system is described and results of its applica- 
tion are presented. 


INTRODUCTION 

Placing a permanently manned space station in low earth orbit has 
been identified as the next major goal of the United States civilian 
space program.^ This station will serve as a multifunctional base for 
scientific and commercial advances in space. It will contain 
laboratories for research in such areas as communication, solar system 
development, material processing, and astrophysics. The station will 
also serve as a platform for satellite repair, thus expanding the life 
span of these expensive space assets and reducing repair costs. 

In-orbit satellite equipment updating will also be possible, thus 
assuring that technological developments are quickly incorporated. 
Additionally, the Space Station will serve as a base for the assembly of 
other space structures which are too large to be first assembled on 
earth and then placed into orbit by the Space Shuttle. 

The National Aeronautics and Space Administration (NASA) has been 
assigned the task of defining a reference configuration for the Space 
Station. From a field of five candidate configurations, NASA has 
selected the "Power Tower" arrangement (see Fig. 1) as the reference 
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configuration. 2 One of the reasons cited by NASA for selecting the 
Power Tower configuration is its extraordinary expansion capabilities. 
This configuration will consist of pressurized modules for habitation 
and work areas, solar arrays for power acquisition, radiators for heat 
dispersion, a docking port for the Space Shuttle, and a light-weight 
truss constructed of carbon epoxy to which all the above are attached. 
Due to the large size and the expected growth of the Space Station, it 
is not reasonable to assume that the structure can be analyzed as a 
rigid body. 

The transfer of crew, supplies, and equipment to the Space Station 
will require frequent docking of a Space Shuttle orbiter with the 
station. It is therefore important that an understanding of the effects 
of docking on. the motion of the Space Station/Space Shuttle system be 
developed. Estimates of these effects on the motion of the proposed 
Space Station configuration are needed to adequately design its attitude 
and translational control systems. Careful investigation of the docking 
process should result in improvements in the reference configuration. 

It is expected that the closing rate of an orbiter with the Space 
Station will be small (on the order of 1.0 ft/sec). However, due to the 
high degree of flexibility of the station, the docking of the orbiter 
may still produce significant deflections of parts of the station. The 
rare possibility of a Space Station control systems malfunction requires 
that docking of the orbiter with an uncontrolled Space Station be 
considered. Furthermore, the dynamic response of an uncontrolled Space 
Station/Orbiter system during docking is of considerable importance from 
the standpoint of control system design. 

Early studies on docking involved investigators such as Williams, 3 
Grub in, Chiarappa, 5 Bray ton, 6 and Cochran and Henderson. 7 With the 
exception of the work done by Cochran and Henderson, these early studies 
were not concerned with the effects of flexibility on docking. In 
considering the effects of flexibility, Cochran and Henderson analyzed a 
system consisting of a rigid target vehicle to which two point masses 
are connected by massless, flexible, extensible rods. A rigid 
rendezvous vehicle was allowed to dock with the target vehicle and the 
effects of the flexibility of the appendages were then analyzed. 

Recently, the problem of spacecraft flexibility has received 
renewed attention. In particular, Levinson and Kane, 8 * 9 * 110 have 
considered the planar docking dynamics of bodies consisting of flexible 
and rigid components. Some of the work done by Levinson and Kane 
involves the docking dynamics of (1) a spacecraft modeled as a 
cantilever beam and a rigid rendezvous vehicle, and (2) a spacecraft 
modeled as a free-free beam and a rigid rendezvous vehicle. 9 * 10 In 
these analyses, the deformation of the structure was represented by 
unconstrained mode shapes which were obtained using a finite element 
approach. It was shown by Hablani 11 that unconstrained mode shapes 
portray the deformation of the structure more accurately than 
constrained mode shapes. Here, "constrained mode shapes" refers to mode 
shapes obtained when one end of the structure is constrained and the 


135 



structure is then caused to vibrate at its natural frequencies. On the 
other hand, unconstrained mode shapes are obtained when the structure is 
completely free. Many other investigators such as Likins, 12 Hughes and 
Skelton 13 and Ho and Herber 1 ** have developed simulation routines and 
modal truncation methods to analyze the effects of flexibility on large 
spacecraft. 

The docking problem considered in this paper differs from the work 
done by Levinson and Kane. Herein, the motion is three-dimensional 
rather than planar and a more complex model of one of the bodies is 
developed and utilized. Unconstrained mode shapes are used to define 
the motion of the structure. These mode shapes were obtained from the 
MacNeal-Schwendler (MSG) Nastran 1 ^ finite-element modal analysis 
routine. They were incorporated in a computer program developed to 
simulate the dynamics of the Space Station before and after docking. A 
particular docking mechanism is not considered; instead, the docking is 
modeled as an impulsive interaction between the Space Station and an 
Orbiter . 


In the following sections the governing equations of the Space 
Station/Orbiter system are derived and simulation results are presented. 
The motions of the Space Station and Space Station plus Orbiter system 
are considered first. Next, the equations governing the impulsive 
interaction between the Station and an Orbiter are derived. Use of the 
complete set of equations to simulate the motion of the Space Station 
before and after docking with an Orbiter occurs is then discussed. 
Representative simulation results are presented and interpreted. 

Finally, conclusions are stated along with suggestions for additional 
research. 


EQUATIONS OF MOTION 


The station is modeled as a hybrid of flexible and rigid components 
in which the base section (lower section containing the modules and the 
docking port) is modeled as a system of rigid bodies and the remaining 
structure (upper and lower keels, booms, and solar arrays) is modeled as 
a collection of flexible bodies. The center of mass of the undeformed 
station is assumed to move in a circular orbit about the earth. Because 
impulsive docking is assumed, the presence in the system of the Orbiter 
after docking is modeled by adding a rigid body to the base section of 
the original Station model. 


As shown in Fig. 2, 
mass element, P, of mass 


the inertial position vector R^ of a generic 
dm, can be written as 


R 3 R + r , 
-p -s - * 


( 1 ) 


where jr includes vectors for both the deformed and undeformed structure. 
If the structure is assumed to be flexible, the vector jr can be 
expressed as 


r 


r + u . 
-u 


( 2 ) 


As shown in Fig. 3, r u is the vector locating the point P in the 

SX 3 Y 3 Z 3 system when the structure is undeformed; and _u defines the 
deformation of the structure in the SX Y Z system. The SX Y Z 

system is tied to the station in such a way that S coincides with the 

station's mass center when the structure is undeformed. The nominal 

orientation of this coordinate system is such that the X 3 -axis is 

tangent to the orbit, the Y 3 -axis is parallel to the boom and points 

towards starboard, and the Z 3 -axis is parallel to the keel and is 

directed radially towards the center of the earth. When docking occurs, 

the SX Y Z system rotates with the Space Station. If the deformation 
s s s 

is defined in terms of unconstrained mode shape vectors, then _u can 

be written as 


i* * l JfcV (3) 

k-1 * 

where ^ is the mode shape of the k th mode, and q^ is the generalized 

coordinate associated with the k c ^ mode. Substitution of Eq. (2) into 
Eq. (1) results in 

R * R ♦ r + u. (4) 

-p -8 -u 


The translational equations of motion are obtained from 



f a dm 



/ dm » 

M p 
s 


(5) 


where f is Che force per unit mass acting on the differential element of 
mass, dm, and M s is the total station mass. 


From Eq. 


(4), R can be expressed as 
~P 


.. »• , o oo 

R * R * to x (r +u) + to x[o) x (r +u)] + 2u> x u + u , (6) 

-p -8 -s -u - -s -s -u - -s - - 

where u) is the angular velocity of the SX Y Z system and a ,,0,, over 

a vector, denotes time differentiation of that vector's X 3 -, Y a - and 
Zg-components only . Substituting Eqs. (3) and (6) into Eq. (5), the 
resulting equation can be evaluated to obtain the following matrix form 
of the translational equations of motion: 
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( 7 ) 


,ex 


N 


M 


R 

-8 


S . I l( »8 + *.V V»k 


k-1 


+ 2 Wk + \\ ] 

where is Che modal momentum coefficient ^ defined as 

“k * / 4^ d ® » 

M 

8 


( 8 ) 


and F ex is the total external force. The direction cosine matrix, g T , 

is also introduced to accommodate the different coordinate systems used 
in writing Eq. (6). Here, C? transforms vector components in SX 8 Y 8 Z 8 

system to corresponding components in the EXYZ system. Thus, Eq. (7) 
represents the translational equations of motion written in the EXYZ 
system. 

The rotational equations of motion about the station's mass center 
are obtained from 

/ j dm * / r x a dm, (9) 

M M ~ ~ P 

s s 

where t is the torque per unit mass about S due to the force per unit 
mass _f. The left side of Eq. (12) is the total external torque about 
the station's mass center; that is, 

/ x dm - T ex . (10) 

M 

s 

By substituting Eqs. (2), (6), and (10) into Eq. (9), and evaluating the 
integrals which appear in the resulting equation, one may rewrite Eq. 

(9) in the matrix form, 

T ex - I l VkSl, ♦ is. ♦ S,I», 

k»l 

* J *lA , k * 2 MA ♦ ( *k * J, 

where I is the inertia dyadic of the station, 
coefficient , H of Eq. (11) is defined as 


*kjVV • (1I) 

The flexibility 
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4 - l 44 *" 

1 s 

and 

4j - I ij4 d " 

M 


Again, the direction cosine matrix, C, has been used to transform the 

components of vectors to a common reference coordinate system (i.e., 
SX s Y s Z s system in this case). 

Equations (7) and (11) constitute six of the six plus N (6+N) 
equations needed to define the motion. The remaining N equations are 
obtained from the equations 

/ 4 i • f dm = / 4 • a dm , j-l,2,...,N, (] 

M J M J P 

s s 

where a p is given by Eq. ( 6 ). Generalized forces, Q j , j-1 ,2,3, . . . ,N, 
may be defined by 

Q: * / 4; • f dm a • f eX dm + / i- • f in dm , (1 

J M J M J M J " 

s s s 

where f ex and f*- n are the "external" and "internal" parts of f. Thus, 
Eq. (13) can be rewritten as 




O 

id x (r +u) + <b. 
-s -u - - r j 


CD X {(D x(r +u) J 
— s l— s -u — 1 


* 2 ij 


e 

0) X u + <fc. 
-s - *j 


u] dm 


(1 


The integrals on the right-hand side of Eq. (15) may be evaluated to 
obtain the following matrix form for Eq. (15): 


T •• 
m.£R 
-j*-s 


+ ( 




N 

I 

j=l 


*kj ■ 


T 

a) R. a) 
-s^j-s 


N 


J * 2 


N T 

f (a $. 




s-^jk^k 


N 

l 

Vs 1 


m u,<i 


jk H k* 


(1 



where 


R. 


/ 6. r 

M J " U 


dm 


(17a) 



u«. 

Jk 



dm , 


dm . 


(17b) 


(17c) 


Since the elastic deformations are assumed small (i.e., less than 
10% of the structural length), then all terms second-order and 
higher-order in q^ are dropped from Eqs. (7), (11), and (16). Similar 

action is taken regarding the angular velocity, u>_ , of the station. 

Also, all products of q^ and w g are neglected. 

By dropping the appropriate terms from Eqs. (7), (11), and (16), 

one may rewrite these equations as follows: 

f x - “A ♦ S T ! I- V k s. * v k i • (18) 

k*l 


f* ■ I iiuSkSi, * i“ s ♦ Vk 1 ' 

k»l 


ip 

Q. * m.CR 
J -J— • 


T. 

Vs 


k a l 


jk^k 


(19) 


( 20 ) 


Explicit expressions for the forces and torques expressed in Eqs. 
(18), (19), and (20) are now developed. 

Forces and Torques 

A satellite in a low earth orbit (nominally 300 km altitude) is 
affected by several forces. For example, the important external forces 
and torques may include those due to gravity, the atmosphere, thrusters, 
other control devices, and solar radiation. In this paper, gravita- 
tional forces and gravity-gradient torques are the only external forces 
and torques considered. 

The gravitational force acting on a differential element of mass, 
dm, is expressed as 
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( 21 ) 


- 8 * ' rT - P ’ 

p 

where u is the gravitational constant of the earth and R is the vector 

o -p 

shown in Fig. 4. 

By recalling Eq. (1), one can approximate R~ 3 by (Ref. 17) 

R ' 3 * jr (I " -|z E * 5 S + H.O.T.] (22) 

F s s 

By substituting Eqs* (l) and (22) into Eq* (21), and then integrating, 
one may show that 


F ex - F « 


P. M _ N 

rJ m A ' T* S I 5k«k 

s s k«l 


3y N 

* “sf SkV*. 


c t t 
2F* ‘ 

s 


: r <P R s - R s C T I + H.O.' 

S 


Next, the contribution by terms containing R~ 5 may be neglected and 

s 

the external force expressed (in the EXYZ system) as 
ex u o T N 

S m ~ rJ t M s 5 s * S I • (24) 

a 88 k-1 * k 

It was noted by Kaplan* 8 that gravity-gradient torques provide 
excellent directional stability for spacecraft in eccentric orbits. 

Thus, in considering a circular orbit (eccentricity equals zero), it is 
reasonable to assume that gravity— gradient torques may be used to 
provide stabilization for the Space Station. The gravity-gradient 
torque can be obtained by crossing the vector jr with the external forces 
defined in Eq. (21) (Ref. 17, pp. 112-119) and is presented here as 


E “ - Vn * h S 9 iss a l- 


* 5? SS.J 

s 


Equation (25) represents a set of differential equations, written 
in the SX s Y s Z s system, that governs the rotational motion of the Space 


Station. 
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The generalized force, Q j , associated with the j c ^ mode of 
vibration can be decomposed into two parts. One part contains the 
contribution due to internal forces, and the other contains the 
contribution due to external forces. Each contribution can be 
formulated as shown in the following set of equations. 

Q j " Q j° + Q j X (26a) 



n-i " ik ’ k * Ji CjkV Ji k j k<lk 


(26b) 




N 


* I 

k*l 


jk q k 


(26c) 


In Eqs. (26), c;fc and kj^ represent the jk th element of the structural 
damping and stiffness matrices. The quantity Yj in Eq. (26) is defined 
as follows: 



(27) 


The complete linearized equations of motion are obtained by 
combining Eqs. (18), (19), and (20) with Eqs. (24), (25), and (26). 

From these expressions, one may show that the equations that govern the 
translation motion may be expressed as 



' "s 5 . * ^Jj 1 ' =k V. * 5k V • <28) 

whereas, the rotational motion is governed by 



■ J/VkSS, * Is. * Vk 1 • < 2,) 

The generalized coordinates are governed by the following 
expression: 



(30) 



T T 

* V 

N 

n 

R C m. 
-s- -j 

+ i m ik q k 
k-1 JK K 

T ~ 

m.CR + 

T • 

U . b) 

N 

+ y m.. 

-j— s 

-J -s 

k-i * 


N 


The 6+N linearized differential equations expressed in Eqs. (28), 
(29), and (30) are sufficient to simulate the motion of the station 
before and after docking. To simulate motion after docking, M s and I 

must be adjusted to account for the presence of the orbiter that is now 
attached to the station. 

DOCKING 

As stated previously, a particular docking mechanism is not 
considered in this paper. Instead, the docking of the orbiter with the 
Space Station is assumed to produce impulsive changes in the kinematical 
variables. In what follows, the docking of the orbiter with the Space 
Station is modeled as rigid body docking with a cantilever beam/rigid 
body assembly. A similar problem was addressed by Levinson and Kane in 
Ref. 10 where they analyzed the planar case of a rigid body docking with 
a free beam. 

The system to be analyzed (see Pig. 5) consists of a single rigid 
body (orbiter) and a cantilever beam/rigid body assembly (Station). The 
system of Fig. 5 can be subdivided into the systems of Fig. 6. 

The Space Station (bodies A and B) is assumed to have linear and 
angular velocities, u s and w s , before docking occurs. After docking, 

these velocities become v and g , respectively. The orbiter' s linear 

velocity will change from u o to v q due to docking. Also, the angular 

velocity of the orbiter will change from « to g . As stated above, the 

orbiter and the Space Station are assumed to couple rigidly during 
docking, and then to rotate together. 

Equations for Docking 

The law of conservation of linear momentum for body A can be 
expressed as 

t, 

/ <Y a - H. )<lm * / </* . (31) 

"* ‘l 

where m. is the mass of body A. The quantities u and v are the 
A “a —a 

"before" and "after" velocities of a differential element of mass, dm, 
and are defined here as 
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Fig. 5. Schematic of Orbiter/Space Station Assembly 



Fig. 6(a). Free-Body Diagram for Flexible Section of 
the Space Station (Body A). 



Fig. 6(b). Free-Body Diagram for Rigid Section of 
the Space Station (Body B) . 



Fig. 6(c). Free-Body Diagram for the Orbiter (Body 0) . 
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S« ■ H, " V 5* * S. x i* ♦ , l Vk 


2, * 8. * 5 a * 9, * £* * I 4^* 

k 3 *! 


• «. 

where and are the "before" and "after" time rate of change of the 

generalized coordinates. Equation (32) may be substituted into Eq. 
(31), and the result integrated to obtain the following form of the 
statement of the law of conservation of linear momentum for body A: 

M , + 

“a - 2.> ♦ "a <3, - V * <3* * IcA> * Jj £k ( ’k • V 


- - f 1 f* dt . (2 

C 1 

The quantity shown in Eq. (33) represents the momentum flexibility 
coefficient * * for body A and is defined as 


4 / A <*" 


The law of conservation of linear momentum may also be applied to 
bodies B and 0. For body B, one has 


mgCv^ * a / (F + f*)dt. 


and for body 0, 


m(v - u ) ■ - / F dt 

o -o -o J 

C 1 


Equations (33), (34), and (35) may be combined to give 


■a ( 2. - »•> * m A < 3s” 2>.) * <3 a * Jca> * ,Mk (q k • V 

k=l 


+ m_(v - u ) + m ( v - u ) = 0 


(36 


The dot product of with Eqs. (31) and (32) combined provides an 

o+ o- 

expression for - q^. 


/ t(v s - u s ) • 4 . + (g s - 0 > s ) x (R a + r*) « 4 . 

TTI J J 

N o+ o _ t. 

+ I 4; * 4i, ( <L " dm =* - / 2 4- • f* dt 

k»l J K K K J 


( 37 ) 


The right-hand side of Eq. (37) is zero because at the point at which f* 
is applied, the deformation is zero. Thus, Eq. ( 37 ) can be rewritten Tn 
matrix form as 


pT(v - u ) - pT R (Q - o) ) + gT (8 - u 

-j -s -s -J -A -s -s -j -s -s 


N 


- 1 - t> 


k -1 


jk"Mc 


(38) 


where 


G. * / r* A. 

'J l ~ 


(39a) 


and 


T 

m * . m f 6 , 4 dm 

jk j 


(39b) 


To account for the N modes of vibration, Eq. (38) can be rewritten in 
the following form: 


E T <! 8 - “ 8 > - E T V5 s - 8.) ♦ S T %- 8 S > - S3 . 


(40) 


T T 

The Nx3 matrices, P and G , are defined as: 


[pT] 


T 

g: 

-1 

• 

-1 

F* 

, g t - 

g! 

-2 

• 

* 3 

-2 

• 

p t 

-N 


g t 

-N 


(41a) 
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(41b) 
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— 
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m 21 
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9 * 
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The law of conservation of angular momentum yields three 
additional expressions. These expressions are obtained from the 
following equations by applying the law of conservation of angular 
momentum to bodies A and B combined (Space Station), and to body 0 
(Orbiter). Equation (42) is obtained when the law of conservation is 
applied to the Space Station. 


<i. * IV 


• (Q - (i) ) + 

-s -8 


N 

■v's* * W * I Sk’k 1 


<2, -a.) 


N N t 

♦ R a * l £ k q k + l S k q k - J 2 [I + Rj. x F]dt 

A k-1 k-1 tj 


(42) 


A similar expression is obtained for body 0 (Orbiter) . 


I 

■o 


• (Q - u ) 
-8 *“8 


- / Z (T + R. x F] 


dt 


(43) 


If one combines Eqs. (35), (42) and (43), then the statement of 
conservation of angular momentum can be expressed in the following 
matrix form: 


<1, ♦ « * fc) 8, ♦ t-A * Jca> * £al(v, - »,) * [j A | * gl i 

* *o ( Sd ‘ V^o ‘ V " ( I, * W s, * io So (44) 

Equations (36), (40), and (44) represent six plus H (6+N) equations 
in nine plus N (9+N) unknowns. The relationship between v g and v q 

accounts for the remaining three equations. This relationship is 
written here as 

^O ’ ?S - ( 5o - V s 8 - (45) 

It can be shorn that Eqs. (36), (40), (44), and (45), when combined, 
yield the following equations that govern the impulsive interaction 
between the Station and the Orbiter. 
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tssr 1 <s T - s t £ a > - « A <i A ♦ i CA > - <So - yi s s 

+ [(tn. + m + m ) E + PM-lp T ] v 

- ((m A * m,) E ♦ E!f 1 E T 1 2 , * “o a 0 
♦ [£ST‘ (§ T - E T £ a > - » A (R a + t CA )] u s ( 46 ) 

u 8 ♦«♦£„-»„ y-y^-y 

* <S A £ * g> if 1 <§ T - £ T £ A >) 3 s 

* ( "a ( 5a * W *!a * "o <?D * V ♦ <5 a £ ♦ S)ST‘£ T 1 I, 

■ [ "a ( 5a * Ica> * K * ( 5a 5 * 2>«' l £ T ) s 8 * % <£d - i Q >a 0 

* ij s * x : A * <i A £ * g> r 1 <s T - £ T i A )] ? s * i 0 <*7> 

3 * S”‘£ T ( i s - y * sr l [g T - £ T » A Ha s - ■ . ) (48) 


One may obtain from Eqs. (46), (47) and (48) the changes in the 
state variables which occur during docking. Consideration is now given 
to the simulation of the motion of the Space Station before and after 
docking occurs. 

SOLUTION TECHNIQUES 

Computer programs were developed and used in conjunction with the 
Harris H-800 minicomputer at Auburn University to simulate the motion of 
the Space Station. The computer programs consist of three steps, the 
first of which simulates the motion of the station prior to docking. 

The second step encompasses the determination of the changes in the 
kinematical variables due to the docking of an Orbiter with the Space 
Station. Finally, the third step simulates the motion of the Space 
Station/Orbiter system after docking has occurred. 

In the following sections, the details of the numerical simulation 
of the Space Station before and after docking are discussed. 
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Mode Shapes 


To describe Che deformations of the flexible structure, 
unconstrained mode shape vectors (eigenvectors) were obtained from a MSC 
Finite Element program available at Auburn University on the IBM 3033 
computer. These mode shape vectors were obtained using the Givens 
(tridiagonal) method, and were normalized with respect to the mass of 
the Station. 

The truss structure of Fig. 1 was modeled as a structure consisting 
of "equivalent" beam elements (see Fig. 7). The equivalence was 
obtained by modeling a section of the truss structure and subjecting it 
to known forces and torques. With the deformations (both linear and 
angular) obtained from this model, the stiffness of an equivalent beam 
was computed using standard beam theory equations. 

Using the equivalent beam model for the station, two sets of 
eigenvectors were obtained. One set contains the eigenvectors for the 
station without the Orbiter and the other set contains the eigenvectors 
for the Space Stat ion/Orbiter system. To obtain the eigenvectors for 
the Space Station/Orbiter system, the Orbiter was modeled as a rigid 
body attached to the rigid base of the Station. 

Numerical Simulation 


Equations (28), (29), and (30) represent a set of six plus N ( 6 +N) 
coupled, second-order differential equations. In the simulation 
process, these equations are numerically integrated using a fourth-order 
Runge-Kutta scheme. In order to efficiently utilize the Runge-Kutta 
scheme, Eqs. (28), (29), and (30) were manipulated to give the following 
expressions. The expression associated with the generalized coordinates 
is shown here as 

[MjJ + [cl a + [K] 3 3 F (49) 

where 

<»i - - f h 

S 

+ [ sr /xa + ^ T] i~ l Hi + *r 2x} » (50a) 

s “ s 

[C] - - 2;U n ]y 


(50b) 
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(50c) 



M * tjij JJj JI3 ... jJj,] , 



(54b) 
(54c, d) 

(59e,f,g) 
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During Che simulation, Eq. (49) is solved for Che highest 
•* 

derivative, 3 , and then integrated to obtain the time rate of change of 

the generalized coordinates, 3 . During each time increment, the 

calculated values of are substituted into Eq. (51), which is 
integrated numerically to obtain the angular velocity and angular 
displacements. A similar process is applied to Eq. (52) to obtain the 
linear velocity components and linear displacements of the station's 
mass center. 

The process described above is continued until docking occurs. At 
the instant of "docking," the changes in kinematical variables are 
obtained by solving Eqs. (46) and (47) simultaneously for the angular 
velocity, and the linear velocity, v , of the station's mass center. 

These results are then substituted into Eq. (48) to obtain the time rate 
of change of the generalized coordinates. 

Motion after docking is simulated exactly as that before docking, 
with the exception of accounting for the presence of the Orbiter. 

Before any simulation may be accomplished, the initial conditions 

o 

on u> , R , q. , q. and the Orbiter s linear and angular velocity must be 
given or calculated. This process is considered in the next section. 

Initial Conditions 


The Space Station's center of mass is assumed to be initially in a 
500 km circular orbit inclined at 28.5° (see Fig. 8 ). Initially, the 
Station rotates at the mean motion for that orbit. For simplicity, the 
angle of the ascending node, 0 , was assumed to be zero, and simulation 
began when the Station occupies the ascending node position. The 
initial conditions of the generalized coordinate were obtained by 
assuming a state of dynamic equilibrium for the Space Station. Using 
these assumptions, the initial conditions on the Space Station state 
variables were calculated and are shown here in Table 1. 

The results discussed in the following sections were obtained using 
these initial conditions. 

RESULTS AND DISCUSSION 

An algorithm was developed to numerically integrate the equations 
of motion of the Space Station and the Space Station/Orbiter system. In 
this algorithm, a fourth-order Runge-Kutta integration scheme with an 
integration time increment of 0.01 seconds was used. The motion was 
simulated over a real time interval of 200 seconds during which docking 
occurred. A separate algorithm was developed to compute changes that 
occur in the kinematical variables of the Station due to impulsive 
docking with the Orbiter. In both algorithms, the structural 
deformations were calculated using the first eight vibrational modes. A 


155 




I>0. 


TABLE 1 


INITIAL CONDITIONS 


Parameter 

R g (EXYZ system) 

x-component 

y-component 

z-component 

(EXYZ system) 

x-component 

y-component 

z-component 

» (SX Y Z system) 
-s s s s * 

x-component 
y-component 
z-component 
9 j i j*l , 2 , . . , f N 


,2, . . . ,N 


Initial Value 

6878 km 
0 km 
0 km 

0 

7.6127 km/s 
0 km/s 

0 rad/s 
0 rad/s 

1 . 1068x1 0“ ^ rad/ s 
-0.19623, -0.41874 
0.25247, -0.52811 
0.23838, -0.13866 
0.10141, 0.22680 

0 , 0 , 0 , . . . , 0 
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third algorithm was developed to compute the mass properties of the 
Space Station, which are shown in Table 2. 

To investigate the effects of docking on motion of the Space 
Station, three approach orientations of the Orbiter were considered. 
These are as follows: (1) An approach along the X-axis in the positive 

X-direction, (2) an approach at 45° to the X-axis in the XY-plane (see 

Fig. 9), and (3) an approach at 45° to the X-axis in the XZ-plane (see 

Fig. 9). For each approach orientation, the closing rate of the Orbiter 
was varied between 0.5 ft/sec and 1.5 ft/sec. 

A comparison was made of the changes that occur in the magnitude of 
the Station's angular velocity when the Station is modeled as a hybrid 
of flexible and rigid bodies to that of the Station when modeled as a 
single rigid body. Figures 10 through 14 show this comparison for 
various approach orientations. It can be observed from these figures 
that flexibility does affect the angular motion of the Station. The 
primary effect is due to the time required for the flexible bodies to 

respond to the impact of the Orbiter with the rigid body to which they 

are attached. This can be called an "inertia" contribution. Although 
the cases of Y-axis or Z-axis approach are probably not realistic, it is 
observed that with either approach the "inertia" contribution of the 
flexible bodies does not produce significant changes in the angular 
velocity of the Station's mass center. Thus, for these approaches, the 
rigid body model yields the larger changes in angular velocity. 

However, for an approach along the X-axis it is observed that the 
inertia contribution of the flexible bodies does significantly affect 
the angular velocity of the Station. Thus, with an X-axis approach, the 
hybrid model predicts a greater change in angular velocity. 

Figure 15 shows the changes in the magnitude of the angular 
velocity of the hybrid model for approaches (1) along the X-axis, (2) at 
45° to the X-axis in the XY-plane, and (3) at 45° to the X-axis in the 
XZ-plane. It can be observed that the XZ-plane approach produces the 
smallest changes in the Station's angular velocity. 

The simulation results of Figs. 16 through 18 were obtained by 
numerically integrating the equations of motion of the Station and the 
Stat ion/Orbiter system. In Eq. (50b), a proportional damping 
coefficient of 0.01 was assumed. Simulation began when the Station 
occupied the ascending node position (see Fig. 8); fifty seconds later, 
docking of the Orbiter occurs. The motion of the Stat ion/Orbiter system 
after docking is simulated for an additional 150 seconds. The results 
represented in these figures are the total displacements of the tip of 
the upper keel (point A) and the center of mass of an upper outboard 
panel (point B) . The deformations of points A and B after docking are 
observed to be in-phase when the Orbiter approaches along either X-axis 
or at 45 d to the X-axis in the XZ-plane. However, when the the 
Orbiter' s approach is at 45° to the X-axis in the XY-plane, the motion 
after docking is no longer in-phase. This is probably due to the 
combined transverse and longitudinal motion of the panels when the 
approach is in the XY-plane. 


TABLE 2 


STATION MASS PROPERTIES 

WEIGHT 

2.6612xl0 5 lbs 


MOMENTS AND 

PRODUCTS 


OF INERTIA 



IXX 

2.0318x10® 

lb-ft 2 

IYY 

1.8704x10® 

lb-ft 2 

IZZ 

2.5434x10® 

lb-ft 2 

IXY 

a 0 

lb-ft 2 

IXZ 

7.4644x10® 

lb-ft 2 

IYZ 

a 0 

lb-ft 2 
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Zs 


XZ-PLANE 


Fig. 9 Approach Orientations for the 
XI- ^nd XZ -Planes. 
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IQsMujgl [(RAD/SBC) x 10 



CLOSING RATE (FT/SEC) 


Fig. 14 Change in Angular Velocity Magnitude vs. 
Closing Rate for an Approach Along the 
Z-axis. 



Figures 16 through 18 also indicate that the maximum deformation of 
point B is substantially greater than the deformation of point A 
(approximately twice as large). Figures 19 and 20 show the maximum 
deformations of points A and B for the various approach orientations and 
closing rates. An approach along the X-axis produces the greatest 
maximum deformations, whereas an approach in the XY-plane produces the 
smallest maximum deformation. 

CONCLUSIONS AND RECOMMENDATIONS 

The results obtained reveal that flexibility is a significant 
factor in the dynamics of the docking of an Orbiter with the proposed 
Space Station. In particular, the changes in the angular velocity of 
the more rigid part of the Station are greatly affected by flexibility. 

Only the first eight modes of vibration were modeled in this 
analysis. Additional modes of vibrations should be considered in 
further studies. However, increasing the number of vibrational modes 
will result in increased computational requirements. 

More general results could have been obtained if the Station's 
payload was considered. In addition, the motion of the crew may have a 
significant effect on the docking dynamics of the Space Station/Orbiter 
system. 

Finally, a suitable control system must be designed to stabilize 
the rotational motion of the Station/Orbiter system. Realistically, 
this can only be done after, or during, a dynamic analysis of a Station 
model that incorporates payload and crew motion. 
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TOTAL DISPLACEMENT (FT) 



Fig.' 17 Displacement Simulation for an Approach at 45° 
to the X-Axis in the XZ-Plane. (Docking Occurs 
at T“50 sec.) 
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MAXIMUM DEFORMATION (FT) 



Fig. 19 Maximum Deformation of Point B 
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